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We analyze four-dimensional, supersymmetric intersecting D-brane models in a toroidal orientifold back¬ 
ground from a statistical perspective. The distribution and correlation of observables, like gauge groups and 
couplings, are discussed. We focus on models with a Standard Model-like gauge sector, derive frequency 
distributions for their occurrence and analyze the properties of the hidden sector. 
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1 Introduction 

The search for string vacua which could provide the Standard Model gauge group at low energy is an 
important aspect of research in string theory in order to make contact with phenomenology. The models 
which have been considered most in the literature are compactifications of the heterotic string and type 
II orientifold constructions (T). The study of type II flux vacua |2jj has led to the conclusion that there 
might exist a huge number (« Googol 5 ) of possible stable vacua. This prospect shows clearly that new 
methods which go beyond explicit constructions of single models have to be developed to deal with this 
huge number. 

One idea that was put forward in 0 and further developed by several authors, suggests to change 
the analysis from a model-by-model construction to a statistical approach. If one were able to detect 
some general behavior in the statistics, independent of the specific model chosen, this might exclude large 
regimes of the landscape of vacua and also give us some hints where to look for realistic models. So far 
most of the statistical analysis has been done in the closed string sector of the theory while there is only 
little work on the gauge sector (for a list of references see @j). In one particularly interesting study 0 
Dijkstra et. al analyzed a huge ensemble of Gepner models to find possible Standard Model candidates. 

In this article we will describe the results of a survey of the distribution of gauge sector observables 
and their correlations for supersymmetric intersecting D-brane models on a T 6 /(Z 2 x Z 2 ) orientifold 
background In earlier work (7J these distributions were estimated by a saddle point approximation, 
later results have been obtained using a brute force computer analysis (4j. This paper is organized as 
follows. In section [2] we will briefly sum up the basic properties of the models we are considering. In 
section [2 we will describe the algorithm we used to compute the results and in section 0 we will present 
the various results from the computer analysis. Finally we will give some conclusions and an outlook to 
further directions of research in section^ 
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2 Models 

2.1 Type II orientifold models 

The class of models we are considering are supersymmetric type IIA orientifold compactifications. The 
orientifold projection O combined with an antiholomorphic involution a, which we will take to be complex 
conjugation in the following, introduces orientifold 06-planes in the background, which carry tension and 
negative RR-charge. To cancel the induced tadpoles one introduces D 6 -branes, which is also very desirable 
from a phenomenological point of view because they will provide us with a realization of low-energy 
particle physics. Both orientifold planes and D-branes, wrap Lagrangian three-cycles 7 r a in the internal 
manifold, which we take to be special Lagrangian in order to preserve half of the supersymmetry. The 
orientifold images of these cycles will be denoted by n' a . 

The homology group //•( (M, Z) of three cycles in the compact manifold M splits under the action of 
flfj into an even and an odd part such that the only non-vanishing intersections are between odd and even 
cycles. We can therefore choose a symplectic basis (a/, /?/) and expand 7 r 0 , 7 r^ and tto6 as 

b 3 /2 b 3 /2 b 3 /2 

Ka = ( X a a I + Y aPl) > K = ( X a a I ~ Y a P 1 ) . *06 = ~Y^L I a I . (1) 

1=1 1=1 Z 1=1 

Chiral matter arises at the intersection of branes wrapping different three-cycles. Generically we get bi¬ 
fundamental representations (N a ,Nb) for two stacks with N a and Nf, branes. In addition we get matter 
transforming in symmetric or antisymmetric representations of the gauge group for each individual stack. 
The multiplicities of these representations are given by the intersection numbers of the three-cycles, 

63/2 

lab := TTa O TT b = £ (Xfo 1 - XlYj) . (2) 

1=1 


The possible representations are summarized in table[T] 


reps. 

multiplicity 

reps. 

multiplicity 

(N a ,N b ) 

lab 

Sym a 

2^aa' -^a06) 

(N a ,N b ) 

lab' 

Anti a 

2 (-^aa' H - -^a06) 


Table 1 Multiplicities of the chiral spectrum. 


The tadpole cancellation condition for k stacks of N a D-branes is given by 

k k 

Y / N a (n a +7r' a )=4no6 - ^ iV 0 X a 7 = l/. (3) 

a =1 a=l 

The supersymmetry conditions are given by the calibration condition for the cycles and an additional 
constraint to exclude anti-branes. 


9(n 3 )| Wo = o, K(f2 3 )k > 0, (4) 

where D 3 is the holomorphic 3-form. Written in the symplectic basis these equations read 

G/2 63/2 

Y ’ fr ='Y^X I a u I >0, where fi~ ^ 3 , «/== / ^ 3 - (5) 

/=1 1=1 ^& 1 J OLi 


In addition to the tadpole and supersymmetry conditions there exists a further constraint from K-theory 
(8l . It is a consistency condition to assure that the Z 2 valued K-theory charge of an Sp(2) probe brane is 
conserved. If this constraint is violated, the model suffers from a global gauge anomaly. 
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2.2 T 6 /(Z 2 x Z 2 ) 

In the sequel we will focus on a specific orientifold construction, a compactification on T 6 /(Z 2 x Z 2 ). We 
will consider only a special class of three-cycles on this torus, namely ’’factorizable” cycles, which can be 
described by their wrapping numbers (m, mj) along the basic one-cycles 7r 2? ;-i , 7r 2l of the three two-tori 
T 6 = lliWj Y’ 2 . To preserve the symmetry generated by Ha, only two different shapes of tori are possible, 
which can be parametrized by 6^ S {0,1 /2} and transform as 

Ha : / 7F2i - 1 ^ n2i ~ l ~ 2bl7T2i . (6) 

( 7t2 i —> —7t2i 

For convenience we will work with the combination 7f 2 j_i = 7r 2 j_i — ^,7r 2 , and modified wrapping num¬ 
bers rhi = rrii + birii. Furthermore we introduce a rescaling factor c := (n? =1 (l — 6j)) to get integer¬ 
valued coefficients. These are explicitly given by (i, j, k € {1,2, 3} cyclic) 

A' 0 = cnin 2 ri3, X 1 = —cnifhjfhk-, Y° = cmim 2 m3, Y l = —crhiTijnk- (7) 


Using these conventions the intersection numbers can be written as 

I ab = (x a Y b - X 6 y a ) (8) 

and the tadpole cancellation and supersymmetry conditions read 


k 3 ^ 

Y / N a x a =L, yy— = o, y, xI ui>^ 

o=l 7—0 Ul 7=0 


(9) 


where we used that the value of the physical orientifold charge is 8 in our conventions and we defined the 
vector L := (8c, {8/(1 — The complex structure moduli Uj can be defined in terms of the radii 

(R^\ R\ 2) ) of the three tori as 

U 0 = R^R^R^, Ui = R^R ( 2 j) R^, i,j,k€ {1,2,3} cyclic. (10) 

Finally the K-theory constraints can be expressed as 

k 1 , fc 

J2 N *Y°e 2Z, N aX: e TL, is {1,2,3}. (11) 

a=1 a=1 


3 Methods 

To do a statistical analysis we need our ensemble of vacua to be as large as possible. If feasible we would 
like to determine all solutions of the tadpole equations, satisfying the constraints from supersymmetry 
and K-theory GJ. Unfortunately an analytic solution to the tadpole equations, which are of Diophantine 
type, is not possible. Even worse, it can be shown that the problem of finding solutions falls in the class of 
NP complete problems (9|- This means basically that there is no way to find an algorithm which generates 
solutions in polynomial time. 


3.1 The algorithm 

To develop an efficient algorithm solving the constraining equations, we split the problem in three parts. In 
a first step possible values for the wrapping numbers X 1 and Y 1 are computed, which fulfill the inequality 


3 3 

o < yy x 1 ih < yy l 1 ih 

7=0 7=0 


(12) 
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that can be derived using (|9j. In a second step we apply a fast algorithm to find partitions of natural 
numbers according to the following equation 

k 

Y J S a = C with S a :=J2 N * U i X a and C:=J2 Ll Ui ■ (13) 

o=l I / 

Finally the terms S a of this partition are factorized into values for N a and X l a using the set of X 1 generated 
in step 1. Only configurations that fulfill the K-theory constraints CD} are taken into account. 

Log(# models) 



Fig. 1 Number of solutions found depending on the absolute value of the complex structure variables Ui. 

Note that in this procedure we have to chose a fixed set of complex structures Uj from the very be¬ 
ginning. To find all solutions we have therefore to loop over all possible values for them. This is exactly 
the point where the bad scaling of our algorithm leads to trouble. The time to compute solutions grows 
exponentially with (7/. Within a total computing time of more than six months on a high performance 
computer cluster we have only been able to compute solutions up to an absolute value for (7/ of 12. The 
only point that saves our day is the experimental fact that the number of solutions decreases exponentially 
with the absolute value of the complex structures, as can be seen in figure^ Being interested in a statistical 
statement for the full ensemble of solutions we do not have to worry about the fact that we have missed 
some solutions with higher complex structures, because their influence on the distribution of observables 
is negligible. 

4 Results 

Using the algorithm described in the last section we calculated a total number of ss 1.6 x 10 8 models. 
Using this data we can now try to analyze the observables of these models. The total rank r of the gauge 
group for example is given by r = N a . The resulting distribution is shown in figure |2(a)] Interesting 
is the suppression of odd values for the total rank. This can be explained by the K-theory constraints and 
the observation that the generic value for Y 1 is 0 or 1. Therefore equation suppresses solutions with 
an odd value for r. This suppression from the K-theory constraints is quite strong, the number of solutions 
reduces by a factor of six compared to the situation where these constraints are not enforced. 

Another interesting quantity is the distribution of U{M) gauge groups, shown in figure |2(b)] We find 
that most models carry at least one (7(1) gauge group, corresponding to a single brane, and stacks with a 
higher number of branes become more and more unlikely. This could have been expected because small 
numbers occur with a much higher frequency in the partition and factorization of natural numbers. 

4.1 Standard Model-like configurations 

The special subset of models we are most interested in consists of those which show some properties 
of the Standard Model or, to be more precise, of the MSSM, since we are considering supersymmetric 
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(a) Rank distribution (b) U (M) distribution 

Fig. 2 Frequency distributions of total rank and U ( M ) gauge groups of all models. 


models only. To realize the gauge group of the Standard Model we need generically four stacks of branes 
(denoted by a,b,c,d) with two possible choices for the gauge groups: U{ 3) a x U (2) b x U( l) c x U(l)d or 
U (3)a x Sp( 2)b x U( l) c x U(l)d- To exclude exotic chiral matter from the first two factors we have to 
impose the constraint that Sym a and Sym b , the number of symmetric representations, have to be zero. 
Models with only three stacks can also be realized, but they suffer generically from having non-standard 
Yukawa couplings. Since we are not treating our models in so much detail and are more interested in their 
generic distributions, we will include these three-stack constructions in our analysis. 

Another important ingredient for Standard Model-like configurations is the existence of a massless 
U(l)y hypercharge. This is in general a combination U(l)y = bed} (l)j of all four U(l)s 

such that the condition JT XiNfYi = 0 is fulfilled. More concretely there are three possible ways to 

construct this hypercharge Qy* = \Q a + \Q C + \Qd , Qy !) = ~\Qa - \Qb or Qy ] = -\Q a - \Qb + 
Qd, where choices 2 and 3 are only available for the first choice of gauge groups. In total we have four 
ways to realize the Standard Model with massless hypercharge 1 . 

Log(# models) 



Fig. 3 Number of quark and lepton generations with (red bars on the left) and without (blue bars on the right) enforcing 
a massless (7(1) . 

The first question one would like to ask, after having defined what a ’’Standard Model” is in our setup, 
concerns the frequency of such configurations in the space of all solutions. Put differently: How many 
Standard Models with three generations of quarks and leptons did we find?. The answer to this question 


1 For a complete list of the realization of quarks and leptons in the different setups please see |4j, tables 2 and 3. 
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is zero, even if we relax our constraints and allow for a massive hypercharge (which is rather fishy from 
a phenomenological point of view). The result of the analysis can be seen in figure 0 This is a rather 
strange result, since we know that models with three families have already been constructed in our setup 
(e.g. [ 1(3). A detailed analysis of the models in the literature shows that all models which are known use 
(in our conventions) large values for the complex structure variables (7/ and therefore did not appear in our 
analysis (see section 134*1 . On the other hand we know that the number of models decreases exponentially 
with higher values for the complex structures. Therefore we conclude that Standard Models with three 
generations are highly suppressed in this specific setup, more specifically the suppression factor can be 
estimated to be about 10~ 9 (cf. 0, section 5.1). 


4.2 Hidden sector 

Although we did not find a nice three generation Standard Model in our analysis, it might still be interesting 
to ask questions about the hidden sector of these models. As has been realized in £9 many of the properties 
of our models can be regarded to be independent of each other, which means that the statistical analysis of 
the hidden sector of any model with specific visible gauge group would lead to very similar results. This 
is indeed the case, as can be seen in figure 0] Moreover, comparing this result with the distributions of 
the full set of models (figlire [2} shows that qualitatively the restriction to a specific visible sector does not 
change the distribution of gauge group observables. 



M 


(a) Rank distribution (b) U (M) distribution 

Fig. 4 Frequency distributions of total rank and U ( M ) gauge groups in the hidden sector of MSSM-models (red bars 
on the left) and MSSM models with massive 17(1) (blue bars on the right). 


The frequency distributions of the hidden sector might be compared with a similar analysis for Gepner 
models carried out in 0 To do a direct comparison with their results, we have to restrict the number 
of branes in the hidden sector to a maximum of three. The resulting distribution turns out to be qualita¬ 
tively very similar to the Gepner model results, although our ensemble of models is much smaller and not 
restricted to three generation models. 


4.3 Gauge couplings 

The gauge sector observables considered in the last paragraphs all belong to the topological sector of the 
theory, in the sense that they are defined by the wrapping numbers of the branes and independent of the 
geometric moduli. This does not apply to the gauge couplings, which explicitly do depend on the complex 
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structures, following the derivation in DU, which in our conventions reads 


1 


Mpi an ck 

2V2 M sKa 


1 



W r>W 


TD\ L J TD 
il| it' 


J2 xl ui, 


/=0 


(14) 


where n a = 1 or 2 for an U(N) or Sp(2N) stack respectively. 

If one wants to perform an honest analysis of the coupling constants, one would have to compute their 
values at low energies using the renormalization group equations. We will not do this here, but instead look 
at the distribution of a s /a w at the string scale. A value of one at the string scale does then not necessarily 
mean unification at lower energies, but it could be taken as a hint in this direction. The result is shown 
in figure[5]and it turns out that only 2.75% of all models actually do show gauge unification at the string 
scale. 

Furthermore we would like to analyze the value of the Weinberg angle sin 2 9 = ay/{ay +oi w ) depend¬ 
ing on the ratio a s /a w and check the following relation between the three couplings, which was proposed 
in DU and should hold for a large class of intersecting brane models 


12 11 , 2 3 1 

— = --1-sm 9 =- - -. 

Giy 3 Gig Giuj 2 GZyt jGig 3 


(15) 


The result is shown in figure[6] where we included a red line that represents the relation | Q3) . The fact that 
actually 88% of all models obey this relation is a bit obscured by the plot, because each dot represents a 
class of models and small values for a s /a w are highly preferred (cf. figure^- 


♦models sirC2 theta 




Fig. 5 Frequency distribution of a 3 /a w in Standard Fig. 6 Values of sin 2 6 depending on ots/a-w. Each dot 
Model-like configurations. represents a class of models with these values. 

Comparing these results with the Gepner model analysis of 0 we find- in contrast to the case of hidden 
sector gauge groups (see section l4~2l . no agreement. The fraction of models obeying G3 was found to be 
only about 10% in the Gepner model case. This might be traced back to the observation that in contrast 
to the topological data of gauge groups we are dealing with geometrical aspects here, which are more 
influenced by the fact that we are working in a large radius regime, while the Gepner model analysis is 
done at small radius. 

4.4 Correlations 

An interesting question that we raised in the introduction concerns the correlation of observables. If differ¬ 
ent properties of our models were correlated, independently of the specific visible gauge group, this would 
provide us with some information about the generic behavior of this class of models. As it turns out, this is 
indeed the case. A nice example is the correlation between the total rank of the gauge group and the ’’mean 
chirality” of a model. The mean chirality is a quantity defined as \ := n Yl,,^b^ a ' h ~ lab), where n is 
some normalization. The difference between the two intersection numbers computes the amount of chiral 
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matter in a bifundamental representation between the stacks a and b. \ measures therefore the amount of 
chirality that is present in a specific model. Looking at the distribution of models for specific values of 
rank and chirality we find that their values are not independent. Moreover, this behavior is generic and the 
shape of figure0does not change if we require a different visible sector gauge group. 



Fig. 7 Frequency distribution of models with Standard Model gauge group depending on total rank and chirality. 


5 Conclusions 

We presented some aspects of an explicit statistical analysis of a specific construction in the class of type II 
orientifold models. An interesting result is that models with three generations of Standard Model particles 
are highly suppressed in this setup. We found that the observables in the hidden sector are independent of 
the specific gauge group in the visible sector. Furthermore we found evidence for a correlation between 
observables. Comparing our results to the small radius analysis of Gepner models in 0 we find agreement 
for topological quantities, but disagreement for more geometrical properties like the gauge couplings, as 
could have been expected. 

It would be very interesting to compare these results with statistical data of other constructions, like 
heterotic compactifications. Equally desirable would be the inclusion of fluxes in our statistical treatment 
and the analysis of grand unified gauge groups like SU (5) in the visible sector. 
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